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Abstract 

We argue that the dynamic extended molecular orbital (DEMO) method may be less 
accurate than expected because the motion of the center of mass was not properly 
removed prior to the SCF calculation. Under such conditions the virial theorem is 
a misleading indication of the accuracy of the wavefunction. 

The first step in any quantum-mechanical treatment of atomic and molecular 
systems is the separation of the motion of the center of mass. The nonrel- 
ativistic Hamiltonian operator with only Coulomb interactions between the 
constituent particles for such systems is of the form Ht = T + V , where T 
is the total kinetic-energy operator and V is the sum of all the Coulomb 
interactions between the charged particles. By means of a straightforward 
linear combination of variables one rewrites the kinetic-energy operator as 
T = TcM + Trei, where Tcm and Trei are the operators for the kinetic energies 
of the center of mass and relative motion, respectively. Then one solves the 
Schrodinger equation for the internal Hamiltonian H = T^ei + V [1-3] . 
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It is well known that the eigenfunctions of Ht are not square integrable. 
For this reason, it is at first sight striking that Tachikawa et al [4, 5] carried 
out their dynamic extended molecular orbital (DEMO) method on the total 
Hamiltonian operator Ht- A question therefore arises: how does this omission 
affect the results of the nonadiabatic calculation of molecular properties?. In 
this letter we will try to answer it. 

Suppose that we try to approximate the energy of the system by minimization 
of the variational energy W = (^Ht'^ = {ip\ Ht \^) / {^\ ^) as in the DEMO 
method of Tachikawa et al [4, 5] . If <y9 depends only on translation-invariant 
coordinates then W = W^ei = (h"^ because (Tcm^ = 0. However, if Lp de- 
pends on the coordinates of the particles in the laboratory-fixed set of axes, 
as in the case of the SCF wavefunction used by Tachikawa et al (see, for 
example equations (10) and (7) in references [4] and [5], respectively), then 
W = (Tcm^ + > W^re/- From the variational principle we know that 
Wrei > Eq, where Eq is the exact ground-state energy of the atomic or molec- 
ular system. Therefore, the use of Ht (instead of H) and a laboratory-fixed 
set of axes for the electronic and nuclear coordinates in if will result in an 
even larger estimation of the molecular energy. 

It is well-known that the SCF wavefunction satisfies the virial theorem [5, 6] 
2 <^T^ = — (V^), but in this case we have a wrong relation because (t'^ = 
Tcm) + (Trei) > (TreiY Therefore, under such conditions the virial theorem 
may be a misleading indication of the quality of the wavefunction. 

Table [1] shows the ground-state energies of some diatomic molecules calcu- 
lated with the internal Hamiltonian operator [2,3] and also the corresponding 
DEMO results of Tachikawa and Osamura [5] who did not remove the motion 
of the center of mass. As expected the uncorrelated SCF energies are greater 
than those in which particle correlation is explicitly taken into account [2,3]. In 
addition to it, we also expect the energy difference AW = W'^^ — W^"^ (where 
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TO and KA stand for Tachikawa and Osamura and Kinghorn and Adamowicz, 
respectively) to depend on the expectation value (Tcm'^ that should decrease 
as the molecular mass increases. In fact, the third column of Table [T] shows this 
trend as expected from the fact that (Tcm^ is inversely proportional to the 
total molecular mass. If this argument were correct then AW would exhibit 
an almost linear relation with the inverse of the mass number A. Fig. [T] shows 
that this is in fact the case for the values of the energy difference shown in 
Tabled 

In order to illustrate (and in some way corroborate) the arguments above 
we consider a simple but nontrivial toy example given by the anharmonic 
oscillator 



Ht = TT^ - 7: TT^ + k{Xi - X2 



,4 



2mi dxi 2m2 dx 



In terms of the relative x = X1 — X2 and center-of-mass X = [7711X1+1712X2) /M 
coordinates, where M = mi + 1712, we have 

"""^ = -2MdX-^ - 2^d7^ + (2) 

where m = 77111712/ M is the reduced mass. The first and second terms in the 
right-hand-side of this equation are simple examples of the Tcm and Trei op- 
erators, respectively, mentioned above. This toy model may seem to be rather 
too unrealistic at first sight but if exhibits some of the necessary features. 
First, it is separable into center of mass and relative degrees of freedom. Sec- 
ond, we can apply simple variational functions of coordinates defined in the 
laboratory-fixed set of axes as well as functions of more convenient relative 
variables. Third, we can calculate the eigenvalues of the relative Hamiltonian 
operator quite accurately, which are useful for comparison. 

To simplify the calculation we resort to the dimensionless coordinates qi = 
Xi/L, where L = [/i^/(mi/c)]^/^, and the total dimensionless Hamiltonian op- 
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erator 



where (3 = mi/m2- Analogously, the relative Hamiltonian operator is given by 



where g = gi — g2 is the translation-invariant coordinate. 

We first consider the variational function ip^{a,q) = exp(— ag^), where a is 

a variational parameter, and the total dimensionless Hamiltonian operator 

([3]). Notice that this trial function depends only on the relative coordinate q. 

The calculation is straightforward and we obtain Wr = 3 ■ 6^^'^{P + l)^/'^/8. 

Obviously, the optimized trial function satisfies the virial theorem (t) 
T,,A = 2 (l>)=6V3(/3+ 1)2/3/4. 



In order to simulate an SCF function of the laboratory-fixed coordinates we 
consider (y9„r(a, gi, ^2) = exp(— ag^ — foga). The calculation is also straight- 
forward and we obtain PV^, = 3 ■ Q^/^{y/p + lf/[8{^+lf/^] > Wr. The 
optimized trial function also satisfies the virial theorem (t) = 2 (v'^, but in 
this case (T) > (Trei) as discussed above. 



Fig. [2] shows Wr, Wnr and an accurate numerical calculation of the ground- 
state energy of the dimensionless relative Hamiltonian operator (jlj) for < /9 < 
1. We clearly appreciate the advantage of using a trial wavefunction of internal 
coordinates, or of properly removing the motion of the center of mass. We do 
not claim that the error in the DEMO calculation of molecular energies [4, 5] 
is as large as the one suggested by present anharmonic-oscillator, but this 
simple model shows (at least) two aspects of the problem. First, that the 
energy calculated by trial functions of the laboratory-fixed coordinates may be 
considerably greater than those coming from the use of relative coordinates if 
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we do not remove the motion of the center of mass properly. And, second, that 
the virial theorem is not a rehable indication of the quahty of the wavefunction 
if it is not based on the relative kinetic energy. 

We can carry out another numerical experiment with the toy model. The total 
mass in units of mi is M/ mi = (1 + Fig. [3] shows that AW = Wnr — Wr 
depends almost linearly on /3/(H-/5) (at least for some values of (3) as suggested 
by the argument above about the actual molecular energies. We appreciate 
that the toy model gives us another hint on the difference between the actual 
molecular energies calculated by Kinghorn and Adamowicz [3] and Tachikawa 
and Osamura [5]. 

Summarizing: if we do not properly separate the motion of the center of mass 
in a calculation of atomic or molecular properties we expect inaccurate results 
unless the approximate trial function depends only on internal, translation- 
free coordinates. Otherwise, the effect of the kinetic energy of the center of 
mass will be a too large estimate of the energy. Under such conditions the 
virial theorem will result in a misleading indication of a supposedly accurate 
wavefunction. These arguments apply to the case in which all the particles are 
allowed to move [5] and may not be valid when some heavy particles [4] (or 
all the nuclei [5]) are considered as merely point charges (a sort of clamped 
nucleus approximation) . 
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Table 1 

Nonadiabatic energies of some diatomic molecules 

Ref. W AW 

H2 

KAOO -1.1640250232 0.111654 
TOOO -1.052371 
HD 

KW66 -1.1654555 
KAOO -1.1654718927 0.102116 
TOOO -1.063356 
HT 

KAOO -1.1660020061 0.0987868 
TOOO -1.068382 
D2 

KAOO -1.1671688033 0.0918650 
TOOO -1.074137 
DT 

KAOO -1.1678196334 0.0885406 
TOOO -1.079279 
T2 

KAOO -1.1685356688 0.0844127 
TOOO -1.084123 



7 



0.20 



0.15 



0.10 



0.05 



0.00 




Fig. 1. AW vs. A ^ for the H2 isotopic series shown in Table [TJ 
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Fig. 2. Ground-state energy of the anharmonic oscillator calculated with the vari- 
ational function of the relative (solid line) and laboratory-fixed (dashed line) coor- 
dinates and the accurate numerical results (circles). 
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3. AW vs. /?/(! + /3) for the ground-state of the anharmonic oscillator. 
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